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Double Michelson/Fabry-Perot interferometer for laser- and displacement-noise- free 

gravitational-wave detection 
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In this paper we demonstrate that a double Michelson interferometer with Fabry- Perot cavities in 
its arms is able to perform laser- and displacement-noise-free gravitational- wave (GW) detection if 
certain model assumptions are met. Assuming the input mirrors of a single Michelson/Fabry-Perot 
interferometer can be rigidly attached to beamsplitter on a central platform one can manipulate 
with interferometer's response signals in a way to cancel laser noise and displacement noise of all 
test masses except the cental platform. A pair of symmetrically positioned Michelson/Fabry-Perot 
interferometers with common central platform can be made insusceptible to the later then, thus 
allowing complete laser- and displacement-noise- free interferometry (DFI). It is demonstrated that 
the DFI response to GWs of the proposed interferometer is proportional to /g„/7, where 7 is the 
cavity half-bandwidth, that is the strongest DFI response allowed by general relativity. 

PACS numbers: 04.30.Nk, 04.80.Nn, 07.60.Ly, 95.55.Ym 



I. INTRODUCTION 

It is widely known that the first-generation laser intcr- 
ferometric GW detectors suffer from the great amount of 
noises of various nature. The major limiting factors at 
low and middle frequencies are seismic noise and thermal 
noise which can be referred to the class of displacement 
noise of the test masses. At high frequencies photon 
shot noise is dominant. In the standard quantum lim- 
ited (SQL) [ij second-generation detectors, being under 
preparation, the cause of SQL is the fluctuating force of 
radiation pressure in the laser beam (back-action noise) 
pushing the interferometer mirrors in a random manner. 
Thus standard quantum limitation also aries due to dis- 
placement noise. 

Each method of suppressing or eliminating of displace- 
ment noise proposed up-to-date is only suited for control 
of only one kind of noise. For instance, active antiseis- 
mic isolation will definitely suppress seismic noise but 
is helpless against thermal noise or quantum radiation 
pressure. From the other hand, quantum-non-demolition 
(QND) schemes of measurements 0, H, 0] are able can- 
celing back-action noise but are certainly not suited for 
dealing with seismic or thermal noise. 

However, recently there has been proposed a revolu- 
tionary new method of displacement noise cancelation 
which simultaneously eliminates the information about 
all external fluctuating forces but leaves certain amount 
of information about the gravitational waves. The major 
idea is to construct such an interferometer that would 
respond differently to the motion of the test masses and 
the GWs. Then the proper linear combination of the in- 
terferometer responses will cancel the fluctuations of the 
test masses leaving a non-vanishing information about 
the GWs. One may find at least two different methods 



proposed up-to-date. 

The first one, described in a series of papers by S. 
Kawamura et al. [1, H, 0] j bears on the distributed na- 
ture of the GWs. This can be best explained from the 
viewpoint of a local observer (or the local Lorentz gauge) . 
In such a reference frame interaction of the GW with a 
laser interferometer adds up to two effects [1, |^. The 
first one is the motion of the test masses in the GW 
tidal force-field. In this aspect GWs are indistinguish- 
able from any non-GW forces since both are sensed by 
the light wave only at the moments of reflection from the 
test masses. We will refer to this as the localized nature of 
the forces acting on the test masses. If the linear scale L 
of a GW detector is much smaller than the gravitational 
wavelength Agw (the so-called long-wave approximation) 
then the effect of the GW force-field is of the order of 



/i(L/Ag 



where h is the absolute value of the GW am- 
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plitude. Relative motion of the test masses, separated 
by a distance L, in any force field cannot be sensed by 
one of them faster than L/c, thus resulting in the rise 
of terms of the order of 0[h{L/\^^Y] describing time 
delays which take the light wave to travel between the 
masses. Second, GW directly couples to the light wave 
effectively changing the coordinate velocity of light, thus 
manifesting itself as an effective refraction index. Light 
wave traveling in such a (boundless) medium acquires the 
information about the GW in its phase gradually, there- 
fore it is a distributed effect having the 0[/i(L/Agw)^] 
order in long-wave approximation. Ultimately, from the 
viewpoint of local observer displacement-noise-free inter- 
ferometry implies the cancelation of localized effects (GW 
and non-GW forces) leaving a non- vanishing information 
about the distributed effect (the direct coupling of the 
GW to fight). 

It was pointed in Refs. 0, 0] that in order the GW de- 
tector to be a truly displacement-noise-free interferome- 
ter it should be also free from optical laser noise since the 
latter is indistinguishable from laser displacement noise. 
Their sum is usually called laser phase noise. Cancelation 
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of laser phase noise in interferometric experiments is usu- 
ally achieved by implementing the differential (balanced) 
schemes of measurements: in conventional interferome- 
ters (such as LIGO) it is the Michelson topology and in 
DFIs proposed in Ref. [3] it is the Mach-Zehnder (MZ) 
topology. 

Although DFI detectors that bear on distributed na- 
ture of GWs allow complete elimination of displacement 
noise, the "payment" for such a gain is the significant 
weakening of the GW response at low frequencies. This 
fact directly follows from the mechanism of noise cance- 
lation: together with displacement noise we also cancel 
GW terms of the /i(L/Agw)° and /i(L/Agw)^ orders, leav- 
ing only the 0[/i(L/Agw)f ]-order terms. In addition, as 
demonstrated in Ref. [3], only 3D (space-based) con- 
figurations of DFIs allow the 2nd-order response, while 
the response of 2D (ground-based) detectors [1] is even 
weaker — it is proportional to h{L/Xgv/)^- For the 
ground-based detectors with L ^ 10"^ m and Agw ^ 10® 
m such a dramatic decrease of susceptibility to GWs 
(^ 10~^) levels all the advantages of DFI in compar- 
ison with traditional detectors dominated by displace- 
ment noise. 

The second type of DFI detectors proposed in Ref. 
[To| . however, is free from this drawback since the mech- 
anism of noise cancelation does not utilize the notion 
of GW distributed nature. The basic element of the 
proposed toy model is a single Fabry-Perot (FP) cav- 
ity pumped through both its mirrors with orthogonally 
polarized light waves. Consider first the pump through 
one of the mirrors. Two responses will correspond to 
the input wave — the refiected wave and the transmit- 
ted one. The latter is proportional to the distance be- 
tween the mirrors which may vary in time under the in- 
fluence of displacement noise and the GW. The reflected 
wave, however, is somewhat different: besides the term 
proportional to the variable length of the cavity it also 
includes the component due to the prompt reflection of 
the incident wave from the input mirror. Thus, reflected 
and transmitted waves carry slightly different informa- 
tion about the GW and the motion of the test masses. 
Combining them properly one is able to exclude displace- 
ment noise of one of the cavity mirrors. In a similar 
way one can deal with the second pump through another 
mirror. Ultimately, the proper linear combination of all 
four responses allows elimination of displacement noise of 
both mirrors leaving a non-vanishing GW signal propor- 
tional to h(L / Xg^)^ in long- wave approximation. Unfor- 
tunately, since the mechanism of noise cancelation bears 
on the effect of prompt reflection which is nonresonant, 
the obtained DFI response is not amplifled with the usual 
resonant factor associated with the accumulation of low- 
frequency signal inside the cavity. 

It may seem that the elimination of mirrors displace- 
ment noise at the /i(L/Agw)°-order contradicts the force- 
like behavior of the GWs at low frequencies. However, 
it has been demonstrated in Ref. [l3| that the loss of 
the resonant gain results in the sensitivity limitation by 



displacement noises of auxiliary optical elements such 
as lasers and detectors. In conventional non-DFI de- 
tectors these noises are negligible since they are finesse 
times smaller as compared to displacement noise of the 
mirrors. In a double-pumped FP cavity these residual 
fluctuations mean that the GW-induced displacements 
cannot be measured absolutely but only with respect to 
some reference test masses, in full agreement with the 
relativity principle. The fluctuations of these reference 
masses will ultimately limit the accuracy of coordinate 
measurements. Therefore, in a double-pumped FP cav- 
ity DFI response is not strictly speaking displacement- 
noise-free: we are able to cancel noise of the mirrors 
(the major "headache" in conventional detectors) but 
not the one of auxiliary optics (lasers, detectors, beam- 
splitters, etc). Another drawback of the double-pumped 
FP DFI is the inability to cancel laser noise since the 
single-cavity-scheme is not balanced. Although one may 
propose several FP-based balanced optical setups, for in- 
stance, modifications of conventional Michelson/Fabry- 
Perot topology, in all such schemes DFI response will be 
still limited by residual displacement noise of some test 
mass(es). This is a fundamental restriction of the model 
in Ref. [ifl]. 

In spite of the undoubtful theoretical "beauty" of DFI 
idea proposed by S. Kawamura et al., such a detector, 
if realized in a full-scale experiment, will surely face cer- 
tain technical difflculties as pointed out in Refs. @, 
For instance, it is hardly possible to have ideal 50/50 
beamsplitters or lossless mirrors. Ultimately, one may 
expect 100 1000-fold reduction of displacement noise. 
In other words, although completely displaccmcnt-noise- 
frce in theory, DFIs of the 1st type will suffer from var- 
ious technical limitations in practice that will eventu- 
ally leave a certain level of residual displacement noise. 
From this viewpoint the theoretical (fundamental) limi- 
tations of the 2nd type (incomplete) DFI schemes may 
not seem so dramatic as one could expect. Eventually, 
in both cases an experimentalist would have to deal with 
certain technical challenges and the only thing that mat- 
ters is the total level of noise left, regardless its nature. 
And though suppression of residual displacement noise in 
DFIs of the 2nd type may turn a more complex problem 
than that in DFIs of the 1st type, the former still have 
a signiflcant advantage — strong GW response at low 
frequencies. 

Therefore, it is reasonable to consider such GW de- 
tectors that would have as simple topologies as possible 
simultaneously allowing as complete displacement noise 
cancelation as possible and having the strongest response 
to GWs that general relativity allows. The latter re- 
quirement can be derived straightforwardly. According 
to relativity principle absolute coordinate and velocity 
measurements are prohibited. This is another formula- 
tion of indistinguishability between the GW and non-GW 
forces at the 0th and 1st orders of L/Agw. The difference 
rises only in the 0[/i(L/Agw)^] order at low frequencies, 
as described above. If L = cr and Agw = c//gw then 
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such a displacement-noise-free GW response will be pro- 
portional to (/gwr)^/i in spectral domain. In addition, 
this response can be further amplified with a FP cav- 
ity, for instance. Then one should expect the rise of the 
(7t)~'^ 3> 1 resonant multiplier, where 7 is the cavity 
half-bandwidth. Ultimately, the strongest DFI GW re- 
sponse allowed by the first principles should be propor- 
tional to ifgw/j)ifgwT)h in spectral domain. 

The optical setup satisfying the requirements of prac- 
tical reasonableness and maximum completeness of dis- 
placement noise elimination, also restricted by condition 
of the strongest response derived above, does not, how- 
ever, immediately follows from some basic principles. It 
is a matter of search at large, limited by practically rea- 
sonable configurations and assumptions. For instance, 
in this paper introducing several model assumptions we 
propose a pair of symmetrically positioned Michelson in- 
terferometers with Fabry-Perot cavities inserted into each 
arm of both interferometers, as a DFI GW detector with 
a reasonably simple optical setup and the strongest pos- 
sible response. 

First, consider a conventional LIGO topology (with- 
out power- and signal-recycling mirrors). Let the end 
mirrors be partially transmittible. In this case an in- 
terferometer will produce three response signals: the re- 
flected (laser-noise-free) one in conventional dark park 
and the transmitted ones in the arms. It is worth noting 
here that certain care is required when calculating the 
responses: each response signal should be evaluated in 
the proper reference frame of the detector that detects 
the corresponding signal Otherwise, unmeasurable 
quantities may arise in the analysis. An experimental- 
ist is able to measure the quadrature components of the 
interferometer responses and record them for further pro- 
cessing. From the set of transmitted signals quadratures 
it is possible, in principle, to construct a laser-noise-free 
linear combination. Therefore, at this stage we may ob- 
tain two signals (quadratures) free from laser phase noise: 
the reflected one and the combined transmitted one. 

Due to the sophisticated frequency dependence of the 
FP cavities responses these two signals can be combined 
in turn to eliminate one of four differential mechani- 
cal degrees of freedoms associated with the test masses 
(beamsplitter, two input mirrors, two end mirrors and 
two end detectors). At this stage we introduce some re- 
strictions into the optical scheme: (i) end detectors and 
end mirrors, and (ii) input mirrors and beamsplitter are 
assumed to be rigidly connected. The practical legit- 
imacy of these assumptions seems questionable and is 
open for criticism, although, no basic principles forbid 
such a gedanken (thought) experiment. Under these re- 
strictions we are left with only two differential degrees of 
freedom, one of which can be eliminated in a combination 
of two laser-noise-free signals. We choose to cancel dis- 
placement noise associated with the differential motion 
of the end platforms. 

Ultimately, due to the symmetry of plane GW wave- 
front we are able to cancel the fluctuations of the cen- 



tral platform (with beamsplitter and input mirrors) 
if the similar interferometer is positioned symmetri- 
cally (see Fig. [3] in the text) and both interferom- 
eters have common central platform (this is another 
gedanken-experiment-supposion). In this case the single- 
interferometer partial DFI responses will have GW term 
of the same sign but the fluctuations of the central plat- 
form will enter with different signs. Adding two single- 
interferometer responses we cancel displacement noise of 
the latter. Then the obtained laser- and displacement- 
noise-free GW response signal turns out to be propor- 
tional to (/gw'T')^ amplified with the cavity resonant gain 



II. SPACE-TIME OF ACCELERATED 
OBSERVER IN THE GRAVITATIONAL- WAVE 
FIELD 

Let us first remind briefiy the "tools" necessary for our 
further considerations. 

As explained in Ref. to obtain physically reason- 
able, i.e. measurable, quantities, calculations should be 
performed in the proper reference frames of the devices 
that produce the corresponding experimentally observed 
quantities. Since they are usually subjected to external 
fluctuative forces, they commit random motions and thus 
we have to deal with their proper non-inertial reference 
frames. Corresponding tools for solving certain bound- 
ary electrodynamical problems in such reference frames 
have been developed in Ref. [llj and utilized in Ref. 
[To| . Therefore, in this paper we will not retell the con- 
tent of these works in detail but will write several useful 
formulas in this section. 

In particular, the space-time of an observer having 
non-geodesic acceleration ^obs(i) along the x-axis and 
falling in the weak, plane, '-I- '-polarized GW h — h{t — 
z/c) propagating along the z-axis takes the following 
form: 



(cdtf 
1 



2 •• 

1 + -^S.ohs{t)x 



y 



h{t~z/c) [cdt-dzf 



(1) 



Conditions of linearized theory |2^obs2;/c^| ^ 1 and \h\ ^ 
1 are assumed to be satisfied for all reasonable x and t. 
Without the loss of generality we may assume y — z ^ 
when considering one-dimensional problems. 

Consider a test mass which in a state of rest (no fluc- 
tuations and no GW) has the coordinate xq with respect 
to the observer (also in a state of rest). If the test mass 
is subjected to some external fluactuative force which 
moves it according to the motion law ^t.m. (t) as seen 
from the laboratory frame (for instance, the one asso- 
ciated with the Earth surface), then its motion law with 
respect to the observer in space-time ([1} is: 



Xt.m.it) ^ Xo + Sxt.niXt), 



(2a) 
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.W-U^W, (2b) 



where ^obs(i) is the observer's law of motion as seen 
from the same laboratory frame. It is assumed here that 
16. m. I, l^bsl ^ \xo\- In fact, Eq. (|2b| is the coordinate 
transformation from laboratory frame to the observer's 
frame. In spectral domain which will be widely used 
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It is also important to take into account the effects 
imposed by the GW and acceleration fields on the elec- 
tromagnetic waves propagating in space-time (H)) . It has 
been derived in Refs. [9|, that the waves propagating 
in the positive and negative directions of the a;-axis can 
be described by the following vector potentials: 



A±(x,t) =A±Q 1 + g±{x,t) + w±{x,t) 
where 



(3) 



a±{x,t) 



9±{x,t) = ifco 



and 



\ ""^^^^ ~c^\ ^^(^^ + i 



h{ti)dh 



w±{x, t) = ifco -S.ohsi't) - ± ^ohs{t) T Cobs(i T x/c) 
c 



Both g±{x,t) and w±{x,t) describe the distributed ef- 
fects: g± is responsible for the direct coupling between 
the GW and the electromagnetic wave and w± describe 
the redshift imposed on the electromagnetic wave by the 
noninertiality of the reference frame. Weak fields a±{x, t) 
describe electromagnetic fluctuations (classical or quan- 
tum). 

In this paper we will also deal with the motions along 
the y-axis. In this case all the formulas remain the same 
but the GW function h{t) should be taken with the op- 
posite sign (this follows from the metric ([1])). 



III. RESPONSES OF A FABRY-PEROT CAVITY 
TO THE GRAVITATIONAL WAVE 

Let us consider the operation of a single Fabry-Perot 
cavity as a GW detector (see Fig. [T]). Cavity is assembled 
of two movable mirrors a and b, both lossless and having 
the amplitude transmission coefficient T, \T\ ^ 1. We 
put distance between the mirrors in the absence of the 
gravitational wave and optical radiation to be equal to 
L. The incident GW h{t — z/c) is assumed to be weak, 
plane, '-I- '-polarized and propagating along the z-axis. 




FIG. 1: Fabry-Perot cavity assembled of two movable mir- 
rors a and b. Cavity is pumped by laser L through mirror a 
with the input wave Ain{x,t) and through mirror b with the 
vacuum-state wave Avac(a;,t). Optical field inside the cavity 
is represented as a sum of the wave A+{x,t), running in the 
positive direction of the x-axis, and the wave A-{x,t), run- 
ning in the opposite direction. The wave refiected from the 
cavity Aaut{x,t) is measured with the (amplitude or balanced 
homodyne) detector Di with the reference oscillation (if nec- 
essary) produced by laser L. Transmitted wave ^Q^^(x,t) is 
measured with the (amplitude or balanced homodyne) detec- 
tor D2 with the reference oscillation (if necessary) produced 
by some local source. 



Then without the loss of generality we assume the cavity 
to be lying in the plane z = along one of the GW 
principal axes, coinciding with the x-axis. 

Cavity is pumped by laser L whose center of mass com- 
mits a fluctuative motion ^l(0 along the x-axis as seen 
from the laboratory frame. Both mirrors a and b have 
associated displacement noise £,a{t) and £,b{t). Finally, 
both the detectors Di and D2 that measure the reflected 
wave and transmitted wave correspondingly fluctuate as 
$Di(i) and CdjC*)- 

To evaluate the response signals of the cavity we should 
perform the calculations in the proper reference frames 
of detectors Di and D2, as pointed above: the reflected 
signal is measured with the first one and the transmitted 
signal is measured with the latter one. 



A. Evaluation of the reflected signal 

Here we will derive the expression for the wave refiected 
from the cavity. Since it is detected by Di we will perform 
the calculation its proper reference frame. This section 
mostly repeats the similar considerations in Ref. p^ . 
therefore we will proceed without detailed comments. 

The origin of the coordinate system is assumed to be 
set up at the center of mass of Di. Therefore, according 
to Eqs. (|2al2bp test masses of the system will have the 
following motion laws with respect to Di: 



2:Di(i) =^Di(t)-eDi(t) = 0, 
XL{t)^^l^{t)-^uAt), 
Xa{t)^Ut)~^DAt), 
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Xbit) ^L+^Lhit) 



-6(t)-eDi(i). 



In last two equations we neglected the small distance 
(compared to the cavity length) between the optical 
bench where laser and detector Di are located and the 
input mirror. 

Let the cavity be pumped by laser L through the input 
mirror a with the input wave 



Ain{x, t) = AinO 1 + g+{x, t) + t) 



X exp< 



ain(a;,t)e-*("°*-'=°^), 



(4) 



Strictly speaking, the argument of xl here should depend 
on a^L itself like t — xi^/c, but since ccl is already the 
quantity of the 1st order of smallncss we can neglect such 
dependence. The vacuum-state pump through mirror b 
can be written as: 



-i uJot+ko{x — L} 



(5) 



Here ain(a;,i) is the "weak" field describing optical laser 
noise of the pump wave and a^acixji) is the "weak" field 
describing vacuum noise in the opposite input port. Re- 
mind, that both the laser and mirror a are located at 
X w 0, where g+{0,t) = ■u;+(0,t) = 0, thus input wave 
does not acquire distributed phase shift when it reaches 
mirror a. 

It is convenient to represent the optical field inside the 
cavity as a sum of two waves, [x, t) and A_ (x, t) , run- 
ning in the opposite directions: 



A±{x,t) =A±o + 9±ix,t) +w±{x,t) 



-i{LOot^kox) 



(6) 



Here a± {x, t) describe the phase shift accumulated by the 
light wave while circulating inside the cavity. 
Output wave reflected from the cavity is: 



outO 



1 + 5-(x, t) + w^{x, t) 
{x,t)e-'^'^'>'+'""'\ 



-i{ujQt~\-kox) 



(7) 



If detector Di is a quadratic amplitude detector then it 
measures the quantity proportional to |v4Q^j(a:, t)p (ne- 
glecting very small terms of the order of fi/cjo). If detec- 
tor is a balanced homodyne detector then it measures the 
quadratures of Al^^{x, t). In this case the reference oscil- 
lation can be produced by laser L. We will call al^^{x, t) 
the reflected signal below. 

To obtain the reflected signal we substitute fields ([H - 
[7]) into the set of boundary conditions (conditions of the 
electric field continuity along the surfaces of the mirrors) 



A+{Xa{t),t) 



{t),t)~ RA^{Xa(t),t) 



A_ixbit),t)= TA^,, {xb{t),t) - RA+{xbit),t), 

{t),t)+TA^Xa(t),t), 

and solve the system with the method of successive ap- 
proximations (see Ref. [I3l)- The required solution of 
the 1st order in spectral domain is: 



TA^o'^iko 



(8) 



The following notations have been introduced: 

'j^2^i{i^0+rt)T 
^ = X - J{2g2i(LUo+n)T ^™0, 

RTe^"^"^ _ R - -Re^'"°^ 

A—O — ^inOj ^outO ^ TP) ^inOj 



7^ = 



I — ^2g2j(a)o+0)T ' 



r2 



J-2 



having the following physical meaning: ^/T^^ describes 
the resonant amplification of the input amplitude Ajno 
inside the cavity, l/'?^g+o describes the frequency- 
dependent resonant amplification of the variation of the 
circulating light wave, Ti. and T are the generalized coef- 
ficients of reflection (from a FP cavity) and transmission 
(through a FP cavity), A_o is the mean amplitude of the 
optical wave inside the cavity running in the negative di- 
rection of the X-axis, ^outo t'^^ mean amplitude of the 
wave reflected from the cavity and is the response to 
GW after a single round trip of light inside the cavity 

It is also very useful to analyze the physical meaning 
of each summand in formula ([H]) : 

1. TZ{a\a — ^inO^fco^L)- This term states that the op- 
tical laser noise Oin is indistinguishable from laser 
displacement noise so both always come together 
and their sum is usually called laser phase noise. In 
spectral domain reflected wave obviously contains 
laser phase noise multiplied by the generalized co- 
efficient of reflection TZ. 

2. Tcvac is the vacuum noise Cvac from the opposite 
input port which is transmitted through the cav- 
ity and comes with the corresponding coefficient of 
transmission T in the refiected signal. 



3. TA_o2iA:o 



/"Cn+o describes the 



(6 + egw)e^""-ea 
part of reflected wave which flows out of the cavity 
and, obviously, represents the phase shift accumu- 
lated by optical wave circulating inside the cavity. 
This phase shift is proportional to the change of 
cavity length due to fiuctuative motion of the mir- 
rors and the change of optical path due to GW 
^gw- Accumulation of the signal is described by the 
resonant factor T^^^j^^- 
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4. Al^^^Q2iko£^a is the phase shift due to the prompt 
reflection from the input mirror. 

5. A^^^^gikQ^D-^ accounts for the motion of detector 
which receives the reflected signal. It is straightfor- 
ward to verify that this summand will be canceled 
in the amplitude detection, evidently meaning that 
the amplitude detector is insusceptible to its own 
motion (since amplitude is measured, but not the 
phase) . 

It is worth noting that the same formula ([8]) one can 
derive without resorting to the non-inertial frame by lin- 
early combining the considerations in (i) the transverse- 
traceless gauge for the GW response and (ii) the labora- 
tory frame for response to the fluctuations of test masses. 
Such a simplification owns to the specific feature of a 
round-trip scheme: the source and the receiver of optical 
radiation are located approximately at the same spacial 
position and thus their clocks could be synchronized al- 
most perfectly. 



B. Evaluation of transmitted signal 

Here we will derive the expression for the transmitted 
signal. Since the source L and the receiver D2 are sepa- 
rated by a large distance, the effect of redshift due to the 
acceleration of detector arises which cannot be obtained 
solely in the laboratory frame. Therefore, we are forced 
to perform the calculations in the proper reference frame 
of detector D2. 

We place the origin of the coordinates at the center of 
mass of D2. Then according to Eqs. (Pa|) and (j2bp other 
test masses move as following: 

a^D.W = CD,W-eD,W = 0, 

XL{t) « -i - ^ Lh{t) + - eD.(i), 

Xa{t) ~-L-\ Lh{t) + Ut) - Cd, (0, 

«6W-eD.(i). 

Again we neglected the small constant distances separat- 
ing laser L and input mirror a, and mirror b and detector 
D2. 

From the viewpoint of detector D2 laser L emits the 
input wave described by the following vector potential: 



Ain{x, t) — A 



inO 



1+ g+{x,t) - g+ -L,t- 



-t- w_|_(a;, t) — —L, t 



x + L 
c 

x + L 



X exp < ^tujQ 



t- 



(9) 



This corresponds to the condition that in the immediate 
vicinity of the source optical wave acquires neither phase 
shift due to the motion of the source, nor phase shift due 
to both distributed effects. 
The vacuum pump is: 



(10) 



It is most convenient to write the wave traveling in 
the positive direction of the x-axis inside the cavity as 
following: 



A+{x,t) = A 



+0 



l+g+{x,t) -5+ y-L,t- 
+ w+(a;, t) — K7+ I —L, t — 



X 



L 



c 

x + L 



c 



X exp-! 



'lUJQ 



+ a+(a;,i)e-4-ot-fco(x+L)J 



(11) 



Such a notation means that the phase is counted from 
X — —L. The motion of the input mirror a will be taken 
into account in the boundary conditions. 

The phase of the optical wave traveling in the negative 
direction will be counted from x = (remind, that both 
g-{x, t) and w^{x, t) vanish at a; = 0): 



A_ [x, t) ^ A_o 1 + 9- {x, t) + w- (x, t) 
+ a„(x,t)e^^(""'*+'^'''^^ 



-i{uot+kQx) 



(12) 



Finally, the transmitted wave with the phase counted 
from x (.g+(0,i) = w+{Q,t) = 0) is: 



^out {x,t)^ ^* l + g+{x,t) + W+ (.T, t) 



t(x,t)e 



— i{LO{]t — kQx) 



^-i{ujQt~kQx) 



(13) 



If detector Di is a quadratic amplitude detector then it 
measures the quantity proportional to |v4*^t(x,i)p (ne- 
glecting very small terms of the order of fi/wo)- If de- 
tector is a balanced homodyne detector then it measures 
the quadratures of Al^^^{x,t). In this case the reference 
oscillation should be produced by some local source. We 
call a*^j(a;,t) the transmitted signal below. 

To obtain the transmitted signal we substitute fields 
([51- fO)) into the set of boundary conditions: 

A+{Xa{t),t) = TAi^{Xa{t),t) - RA_{Xa{t),t), 

A_{xb{t),t) = TA^^{xb{t),t) - RA+{xb{t),t) 
Al^^ {xb{t),t) = RA^,, (xb{t),t) + TA+ {xb{t),t). 

The required 1st order solution of this system in spectral 
domain is: 
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(6+ew)e 



r.t. \ „zf2r 



(14) 



where 



A, 



T 
V: 



^inT)Lh~- -^(1 



A 



inO, 



4* 



X2 



■A 



inO- 



Here is the forward-trip GW response, A+o is the 
mean amphtude of the wave travehng in the positive di- 
rection of the a;-axis inside the cavity and A^^^q is the 
mean amphtude of the transmitted wave. 

Again it is useful to analyze the physical meaning of 
each summand in formula (1141): 



1. T 



Au 



oiko (Cl + 4w - i^^^Da) describes 
laser noise transmitted through the cavity. The 
summand with ^g4' ~ 'i^T^Y)2 rneans that the initial 
phase shift due to laser displacement is additionally 
redshifted with the GW and the motion of observer 
(detector D2), because GW and acceleration fields 
change the rate of laser clock with respect to detec- 
tor D2 clock. It should be also noted that even the 
amplitude detector will be susceptible to — iSlr^Da 
(Cd2(^)''" in time domain), since phase modulation 
is transformed into amplitude modulation in a FP 
cavity. 

2. TZa^s.c is the vacuum noise reflected from the cavity 
into the transmitted port. 



3. B?Te^^'^^''A+o2ih 



describes the total variation of the phase accumu- 
lated inside the cavity. 

4. Al^^^Qiko^u2 accounts for the displacement of the 
receiver. If detector is the amplitude one, this term 
becomes unmeasurable. 



B. 



B 



out 



BS 



BS 

^ out 



!— Hi 

I ^ out » 



A. 



Pi 



FIG. 2: A Michelson/Fabry- Perot optical setup. Two Fabry- 
Perot cavities aibi and 0262 are inserted into Michelson in- 
terferometer horizontal and vertical arms correspondingly. In- 
terferometer is pumped by laser L through beamsplitter BS 
with the input wave Ai^{x,t). Beamsplitter produces to waves 
Ain{x,t) and Bin{y,t) which pump horizontal and vertical 
cavities respectively. Interferometer is tuned to the dark-port 
regime so that both reflected waves Al,^t{x,t) and B^^^{y,t) 
destructively interfere at the beamsplitter and all mean power 
returns towards laser L (not shown in the Fig.). The signal 
part Cout {y, t) containing the accumulated phase shift pene- 
trates into the dark port and is incident on detector Do which 
operates as a balanced homodyne detector with the reference 
oscillation produced by laser L. The dark port of detector 
Dc also produces the vacuum pump Cvac(t/,i). Transmitted 
waves Al,^i{x,t) and B^^t^{y,t) are detected with detectors Da 
and Dfl correspondingly which may operate as amplitude or 
balanced homodyne detectors. In the latter case reference os- 
cillations are produced by some local sources. Both detector 
ports also produce vacuum pumps Avax:{x,t) and _Bvac(j/, i)- 



IV. MICHELSON/FABRY-PEROT TOPOLOGY 

The major disadvantage of a single cavity-based GW 
detector is the significant level of laser noise which dom- 
inates over other noises in practice. To cancel laser noise 
one should implement a balanced optical setup, for in- 
stance, a Michelson interferometer tuned to dark-port 
regime. 

A. Response signals of a Michelson/Fabry-Perot 
interferometer 

Let us consider a Michelson-type interferometer with 
FP cavities in its arms (see Fig. Laser L which ran- 



domly moves along the x-axis as Cl(^) pumps the inter- 
ferometer with the input wave Ai^{x, t). Upon arrival to 
50/50-beamsplitter BS which may fluctuate along the x- 
and j/-axes as Cbs(0 and 77Bs(i) respectively, the input 
wave is splitted into two waves Ain(a:,t) and Bia{y,t) 
which pump horizontal and vertical arms respectively. 
Fabry-Perot cavity in the horizontal arm is assembled of 
two mirrors ai and bi which fluctuate along the x-axis as 
^ai {t) and (t) . The similar cavity in the vertical arm is 
assembled of mirrors 02 and 62 which fluctuate along the 
y-axis as 7702 (i) and rni^{t). Both the cavities may pro- 
duce reflected and transmitted waves. Reflected waves 
Al^^^{x,t) and B^^^{y,t) return towards the beamsplit- 
ter and interfere. Assume the interferometer is tuned to 
dark port regime. This means that the reflected waves in- 
terfere destructively and the mean optical power returns 
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towards the laser. However, the weak time-dependent 
(signal) part Cl^^{y, t) penetrates into the dark port and 
falls on detector Df which fluctuates along the y-axis as 
t7Dc(^)- Lst the latter one operate as balanced homo- 
dyne detector with the reference oscillation produced by 
laser L. Dark port of Dc also produces the vacuum pump 
Cvac{y,t). Transmitted waves A*^t(x,t) and Bl^^{y,t) 
are measured with the corresponding detectors: in hori- 
zontal arm it is and in vertical arm it is Ds moving 
randomly as ^^0^(0 ^^'^ 'nDsi'^) correspondingly. Both 
detectors may operate either as amplitude or homodync 
ones. In the latter case reference oscillations should be 
produced by some local sources (see below). Both de- 
tector ports also produce vacuum pumps Avudx^t) and 

-Bvac(y,0- 

Let us now write the equations for input and output 
waves in an interferometer. We will not write the ex- 
pressions for the waves explicitly; one can obtain them 
straightforwardly making obvious changes in the formu- 
las from the previous section. At beamsplitter the rela- 
tion between the input waves is as following: 

Bin{yBs{t),t) ^ ^AL{xBs{t),t) Cvac (yes (i),t)- 

The relation between the reflected waves is: 

Cout (yes (0 , <) - ^ Blut (yes (i ) , - ^ ^out (a^BS (<),*) ■ 

At this stage we do not define the reference frame, there- 
fore, fields and coordinates of the test masses should be 
specified explicitly for this or that frame. 

These equations can be solved straightforwardly. How- 
ever, we do not need to do this since we already know the 
solution for a single cavity ([5]). First, we need to write 
explicitly the expressions for weak reflected fields a^y^. 
and 6qjjj . To do this we use the first pair of BS boundary 
conditions to obtain in spectral domain: 



1 D _ 1 ^ 



1 



1 



(15a) 
(15b) 



&in = ^ aL + Cvac + Ai,oiko{SxBS - ^Vbs), 

(15c) 

Now let us specify the reference frame. Since both re- 
flected waves will ultimately end up at detector Dc, it 
is necessary to work in its proper frame. Another way 
is to use the laboratory frame which implementation is 
justified with the round-trip situation: laser L and de- 
tector Dc can be approximately considered as located at 
the same spatial position. In any case, substitution Eqs. 
p5a| - [T5c)) into formula ([5]) we obtain: 



1 



V2 



7^(aL 



V 



TA^Q2tko ^ 



^out = '^[«L + Cvac - ^LO«fco(^L " CbS + ??Bs)] + TK 



TB^o2iko 



Here we assumed that the '-|- '-polarized GW is perfectly 
aligned along the interferometer arms. We also neglected 
the term proportional to 77dc since the signal wave does 
not include "strong" mean component. For simplicity let 
us assume that both the cavities have equal detunings 
and bandwidths. This results in A^q = i?_o a-nd A^^jq = 

^outO- 

The boundary condition for reflected waves dictates 
that: 



1 



outO ~ onto 

71 ~ 



1 



V2 



4'' 

^outO' 



71 71 - SyBs)- 



Substituting here the obtained expressions for a^^^ and 



6q^j we obtain: 



TZcv 



V2 



vac 



1 

71 



TA^o 



A'' 

^outO 



+ [^outO + T^^Lo] ifco(^BS - ?7BS). (16) 

One can note that the obtained signal is very similar 
to the one of a single cavity. Namely, the following 
degrees of freedom are equivalent from this viewpoint: 

ia <-> ?7a2 - > 6 ^ "^62 " 61 > ^ ?7bs " Cbs ■ The lat- 
ter relation means that the beamsplitter effectively cuts 
all laser phase noise introducing, however, its own dis- 
placement noise. 

In an experiment homodyne detector Dc measures the 
quadrature components of c^^^. However, keeping this in 
mind, we will deal with the field amplitude itself, since 
calculations with quadratures result in very cumbersome 
formulas, while not changing the physical meaning of the 
ultimate results. 

Now let us derive the transmitted signals. Since they 
are detected by two different devices, each of the signals 
should be calculated in the proper reference frame of the 
corresponding detector. Keeping in mind that fees and 
(52/BS should be explicitly specified for each of these refer- 
ence frame, we substitute Eqs. (|15al ~ [T5c)) into formula 
(fT4l) and obtain: 



1 



V2 



r 
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vac ; 

= ■-i=T(aL + Cvac) + 7^6vac 



Two different regimes of detection of transmitted signals 
are possible. 

1. Resonant regime. Let both cavities be tunes to res- 
onance. In this case amplitude detector measures 
^outo(aout)^ + (^outo)^ flout coinciding with the am- 
plitude quadrature ~ a*^t -I- (flout )^ measured by 
the homodyne detector, since A^^^q is a pure real 
quantity. Resonant regime means that we are tuned 
to the peak of the resonant curve. At this oper- 
ating point variation of optical wave amplitude is 
very weak (~ Q,^). Therefore, neither amplitude 
detector nor homodyne detector measuring ampli- 
tude quadrature can be used. Instead, one should 
use the homodyne detector that measures the phase 



quadrature 



1-)'^. However, in this case 



all the homodyne detectors should be synchronized 
with enough accuracy so that they have equal mean 
phase. Otherwise, different detectors will measure 
slightly different quadrature components. 

2. Non-resonant regime. In this case we are tuned to 
the slope of the resonance curve and variation of 
amplitude of the optical wave is Vl. Amplitude 
detection can be used then. Its major advantage 
is that it does not require synchronization between 
different amplitude detectors. One can also use ho- 
modyne detectors to measure amplitude or phase 
quadratures. It this case, however, synchronization 
between detectors will be required. 



ear combinations between them. Let's consider the co- 
sine quadratures of the signals for definiteness. A simple 
subtraction of 0*^^. quadrature from 6*^^ quadrature, ev- 
idently, cancels the term containing laser phase noise. 
Similar, one can operate with the sine quadratures. This 
can be thought of as a possible method of laser noise can- 
celation from transmitted waves. In the case of reflected 
waves elimination of laser noise takes place at the level 
of interference of field amplitudes and further recording 
of laser-noiseless field in the form of quadratures. In the 
case of transmitted waves we first record the quadratures 
containing laser noise and then linearly combine them 
to produce the laser-noise-free quantity. However, the 
change in a sequence of procedures (to combine first and 
then record or first record and then combine) does not 
introduce any meaningful physical difference, therefore, 
in theory we may operate with transmitted field ampli- 
tudes without the need to perform cumbersome calcula- 
tions with their quadratures. 

Keeping this in mind we construct the following com- 
bination of the fields: 

rft - J- 6* —n^ 

"out "out out 

- ^ TAi^oiko [77BS - ^BS - i^ririBg - £.Ba)] 
v2 ■'ujo+ii 



1 



— i?2Te3'"«M+o 2tko ^^2-2 

v2 ■'ujo+n 



+ ^outO*^o(?7Df 



(17) 



Here we assumed that A+q = B+q due to equal character- 
istics of the cavities. Again, one can establish the equiva- 
lence between this differential signal with the transmitted 
signal (|14p of a single cavity. 



B. 



Cancelation of displacement noise of the end 
mirrors in a single interferometer 



Let detectors and be homodyne detectors for defi- 
niteness so that they measure the quadratures of a\^^ and 
&*ut correspondingly. The reference oscillations should 
be produced by some local sources, for instance, lasers 
that have the same carrier frequency wq and are synchro- 
nized with laser L. As usually required for the homodyne 
detectors, the amplitudes of these local oscillators are 
assumed to be much larger than the mean output ampli- 
tudes A*utQ and Bouto- ^^^^ <^^se we can neglect their 
intrinsic noises (laser noises). However, they are required 
to be synchronized, i.e. have the same homodyne phase 
in order to measure identical quadratures. 

Once the quadratures are measured they can be stored 
in a computer memory and later processed. For in- 
stance, an experimentalist may produce any desired lin- 



Now we have two laser-noise-free signals ([T5)) and p7)) 
which can be combined to cancel one of the fluctuative 
degrees of freedom. Since there are four such quantities, 
Vt^b - Cda: '7&2 - 61, "^02 - Cai and t^bs - Cbs, we should 
somehow suppress two more degrees of freedom by hands 
(the last one will be eliminated by the additional inter- 
ferometer, see below). We introduce the following model 
assumptions (see Fig. [5]): 

1. Both the input mirrors are rigidly attached to 
beamsplitter, i.e. — Cbs and 7703 — riBs- The 
composite mass will be called platform Pes with as- 
sociated fluctuative degree of freedom rjp-g^ — ^p^g . 

2. Detectors and are rigidly attached to the 
end-mirrors bi and 62 respectively, i.e. = 
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and rjh^ = tjub ■ Corresponding platforms will be 
called Pi and P2 and their differential degree of 
freedom rip^ — Cpi- 

The realizability of these requirements remains an open 
practical question. 

Let us now substitute the introduced relations between 
displacements into signals p6)) and p7|) and rewrite them 
in terms of the input amplitude Alq- 



TZcv 



Q-vac ) 



1 i?T^e^*^°^ 
1 i?T^e^*""^ 

- — AwtkoivBS - &s)(l + e'^'^") 



2 T2 T2 

1 



1), 



^out '^Cvac ^ '^(^vac ^vac) 



1 y2giwoT 

:^ ^2 T-2 ^Loifco(?7P2 -Cpi) 

X (1 + i?2g2^(a;o+J^)r ^ jf^^rj^ g*"^) 

1 T^e*'^"^ 

2 7-2 7-2 

1 r^e'"»^ Wo 



2 2za;or \ zOr 



rj^(e^^^^-l) + i?^e^*'^°^(e 



■2 2zcJoT / 2inr 



1) 



From these signals we can exclude either tjbs — Cbs or 
T7P2 — ^pj^ . The following linear combination cancels the 
later quantity: 

1 RT'^e^''^"'' 



vac 



2 7-2 7-2 



^Loi^o(?7BS - Cbs) 



q-2 /I „2ir2T>i 



1 RT^e^''^"'' 



2 ,'1 I g2ir2r -jgiOr 



2 7-2 7-2 



(18) 



where s™'^ is the combined vacuum noise. 

It is straightforward to verify that in the long-wave 
(fir ^ 1) and narrow-band (T^ = 2jt <C 1, where 7 



is the cavity half-bandwidth) approximations signal si 
reduces to: 

5i — 2Cvac V^(&vac ^vac) 



■A 



"f — iS {'J — iS — if2)r 
X ikoinr f ( 77BS - ^BS + ^ i/i ) , 



where d is detuning from resonance. 



C. Cancelation of displacement noise of 
beamsplitter platform in a double interferometer 

The tidal structure of metric ([T]) immediately suggests 
the method of cancelation of beamsplitter platform noise. 
Consider a scheme with two Michelson/Fabry- Perot in- 
terferometers having common central platform (see Fig. 
131) ■ Let us assume that we have eliminated displacement 
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FIG. 3: A double Michelson/Fabry- Perot interferometer hav- 
ing common central platform. 



noise of the end-platforms of the second (left-bottom) 
interferometer and obtained the signal S2 containing the 
fluctuations of central platform and the GW. This signal 
can be evaluated straightforwardly from formula (fT8|l re- 
placing (77BS - Cbs) -{ms - Cbs) and keeping the 
GW function h unchanged due to the symmetry of GW 
wavefront. Ultimately, adding S2 to si we obtain signal s 
free from displacement noise of the cental (beamsplitters) 
platform: 



2in 



11 



2q-2 2iClT 



(19) 



Here s™'^ describes total vacuum noise in both interfer- 
ometers. In long-wave and narrow-band approximations 
we obtain: 



s = 2 



.(1) 



.(2) 



V2 bill 



^vac vac vac 



^ — i5 



■A 



LO 



j:—^,^k,m^\Lh. (20) 



Here vacuum fields with upper index (1) denote the vac- 
uum fluctuations in detector ports of the first interferom- 
eter and vacuum fields with index (2) denote the vacuum 
fluctuations in the corresponding ports of the second in- 
terferometer. 

It is convenient for methodological purposes to com- 
pare the susceptibilities to GWs of the considered in- 
terferometer and of the conventional interferometer with 
Michelson/Fabry-Perot topology (without any recycling 
mirrors) which response is described by the formula 



1 



LO 



X ikoLh{Q) 



sin Qi 



(21) 



Let the readout schemes in both interferometers register 
the following quadrature(s): 



V2i ' 



where 5(57) is given either by formula ((19)) for 
displacement-noise- free double Michelson/Fabry-Perot 
topology (without s™'^ term), or (pij) for conventional 
Michelson/Fabry-Perot topology. To compare the GW 
sensitivities we define the transfer function H{^1) as the 
ratio of GW response quadrature s{fl) to Ai^oikoLh{n). 
We plotted the absolute values of both transfer functions 
in Fig. m for the following parameters (most close to Ad- 
vanced LIGO): L = 4 X 10^ km and R = 0.997. For 
comparison we chose two values of detuning for each sys- 
tem: 6/271 = 0, 100 Hz. 

Fig. d] clearly demonstrates strong GW suppression at 
low frequencies according to the predicted fi^-law in Eq. 
(PO)) . However, at higher frequencies ~ 10** Hz both tra- 
ditional and DFI topologies acquire approximately equal 
level of GW susceptibility. 



V. DISCUSSION 

Let us briefly summarize all the essential model as- 
sumptions that we used for our gedanken experiment. 

1. For a single Michelson/Fabry-Perot interferometer 
we assumed that the end-photodetectors are rigidly 
attached to the end-mirrors and input mirrors are 



■M/FP,8/2jt=0 Hz 
■M/FP, 8/271=100 Hz 
DPI M/FP,8/2jt=0 Hz 
DPI M/PP, 8/271=100 Hz 
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QJ2n, Hz 

FIG. 4: Absolute values of GW transfer functions for conven- 
tional Michelson/Fabry-Perot interferometer {5/2tt = 0, 100 
Hz), and for double Michelson/Fabry-Perot DFI {5/2n = 
0, 100 Hz). 



rigidly attached to the beamsplitter. For a pair of 
interferometers we assumed that both beamsplit- 
ters and all the input mirrors are mounted of the 
common central platform. Although these assump- 
tions do not contradict any fundamental princi- 
ples, the question of their practical realization is 
highly questionable, at least for the ground-based 
facilities. In addition, even such a composite plat- 
form does not allow full cancelation of its inter- 
nal thermal noise (our substitution (r/BS — ^bs) ^ 
~('7bs ~ Cbs) is valid for displacement of the cen- 
ter of mass only). However, in principle, one may 
think of constructing a space-based interferometer 
with arm-cavity lengths of several hundreds meters 
or kilometers which will be most sensible to GWs 
at frequencies below 1 Hz. At such low frequen- 
cies the model of rigid platform may look more at- 
tractive from the experimental point of view than 
that at higher frequencies in the Earth-bound envi- 
ronment. Another way to soften the complexity of 
the optical scheme is to "squeeze" geometrically the 
additional interferometer so that both interferom- 
eters share the same beamsplitter. Nevertheless, 
this does not cancel the requirement that all the 
input mirrors are attached to the beamsplitter. 

2. If the end-detectors of the cavities operate as ho- 
modyne detectors then the local oscillators are as- 
sumed to be present. The frequency of these lo- 
cal oscillators should coincide with the carrier fre- 
quency of the source laser, so they should be kept 
synchronized with it. In addition all the homodyne 
detectors themselves should be synchronized with 
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each other in such a way that they have identi- 
cal homodync phase, otherwise, they will measure 
different quadrature components. The use of am- 
plitude detection looks much more attractive from 
the practical point of view. 

VI. CONCLUSION 

In this paper in a form of gedanken experi- 
ment we have analyzed the operation of a double 
Michclson/Fabry-Pcrot interferometer performing the 
laser- and displacement-noise-free gravitational- wave de- 
tection. It has been demonstrated that if certain model 
requirements met (input mirrors and beamsplitters can 
be rigidly mounted on a single platform and end detec- 
tors can be rigidly attached to the end mirrors) it is pos- 
sible to construct such a linear combination of interfer- 
ometer responses (their quadratures) that produces the 



strongest displacement-noise-free response to the gravi- 
tational wave allowed by general relativity. Namely, the 
DFI response function turns out to be proportional to 
{fgwTY /{it), where /gw is the GW frequency, r = L/c 
with L being the length of interferometer arms and 7 is 
the cavity bandwidth. However, the question of prac- 
tical realizability of our model assumptions is open for 
criticism. 
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